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A stochastic process, when subject to resetting to its initial condition at a constant rate, genetically reaches 
a non-equilibrium steady state. We study analytically how the steady state is approached in time and find an 
unusual relaxation mechanism in these systems. We show that as time progresses, an inner core region around 
the resetting point reaches the steady state, while the region outside the core is still transient. The boundaries of 
the core region grow with time as power laws at late times with new exponents. Alternatively, at a fixed spatial 
point, the system undergoes a dynamical transition from the transient to the steady state at a characteristic space 
dependent timescale t*{x). We calculate analytically in several examples the large deviation function associated 
with this spatio-temporal fluctuation and show that genetically it has a second order discontinuity at a pair of 
critical points characterizing the edges of the inner core. These singularities act as separatrices between typical 
and atypical trajectories. Our results are verified in the numerical simulations of several models, such as simple 
diffusion and fluctuating one-dimensional interfaces. 
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Introduction .— Consider a stochastic process evolving un¬ 
der some given dynamics that does not lead to a time- 
independent stationary state. The prototypical example be¬ 
ing the position of a diffusive particle—which has a Gaussian 
distribution with its width growing as the square root of the 
time, inferring the absence of a steady state. Now, imagine 
that by another mechanism, the dynamics is repeatedly being 
interrupted and recommenced at random times from the ini¬ 
tial condition. A general interesting question is; how does 
such a stochastic resetting affect the temporal evolution of the 
system? 

Examples of stochastic resetting are found in a wide vari¬ 
ety of situations. For example, while looking for a friend in a 
crowded tourist place, after an unsuccessful search over some 
time period, one often returns to the most favorite hangout 
and restarts the search process again. In the ecological con¬ 
text, the animal movements are often modeled by stochastic 
processes EEl. The movements of foraging animals usu¬ 
ally involve local diffusive search for food, interrupted by long 
range non-local resetting moves to relocate to other areas as 
well as to return to their nests, followed by restart of the search 
process 0. For instance, the mobility data of free-ranging 
capuchin monkeys is described quite well by a model of ran¬ 
dom walks with preferential relocations to places visited in 
the past a. Similar notions can be also found in biological 
contexts, where organisms use stochastic resetting or switch¬ 
ing between different phenotypic states to adapt to fluctuating 
environments 00. In computer science, random walks with 
stochastic restarts turns out to be a useful strategy to optimize 
search algorithms in hard combinatorial problems GoHni. 

The stochastic resetting mechanism has been shown to have 
rather rich and dramatic effects on the diffusion process M- 
[m, as well as on long range jump processes such as Fevy 
flights ifTOll . While in the absence of the resetting, the diffu¬ 
sion in free space does not have a stationary state, a non-zero 
rate of resetting to a fixed position leads to a non-equilibrium 
stationary state (NESS) with non-Gaussian fluctuations, in the 
time t —t limit. Similarly, an extended system like a fluc- 
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FIG. 1. (Color online) NESS gets established in a core region around 
the resetting center O whose right and left frontiers (t) grow with 
time. Outside the core region, the system is transient. 


tuating interface, evolving under its own dynamics and is re¬ 
set at a constant rate to its initial configuration, approaches at 
late times a NESS with a nontrivial interfacial height distribu¬ 
tion ll20l . Resetting induced NESS has also been studied in 
other many body systems such as coagulation-diffusion pro¬ 
cesses ED. 

While the mechanism by which the stochastic resetting 
leads to an eventual NESS has been well understood in the 
above studies, the approach to the NESS in such systems is 
yet to be addressed. The goal of this Rapid Communication is 
to study this temporal relaxation to NESS in a wide variety of 
single particle, as well as many body interacting systems such 
as fluctuating interfaces. It is useful first to summarize our 
main results. We compute exactly the time-dependent prob¬ 
ability distribution function (in presence of a resetting rate r) 
Pr{x,t) —where x may represent the position of a single par¬ 
ticle undergoing a stochastic motion with resetting to the ori¬ 
gin, or it may represent the height of a fluctuating interface 
at a fixed point in space with resetting to the initial condi¬ 
tion. As f —oo, the system reaches the NESS characterized 
by P, {x,°°). In systems without the x^ —x symmetry, our so- 
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lution at finite t shows that as time progresses, the NESS gets 
established in an inner core region [—i*_(f),i*+(f)] around the 
point where the system is reset (which we take it to be the 
origin O) —outside this core regime, the system is still tran¬ 
sient (see Fig. [^. The frontiers ^±(f), separating the inner 
NESS regime from the outer transient regime, typically grow 
with time as power laws, thus establishing NESS on larger 
and larger length scales. This phenomenon is conveniently 
captured by probing the PDF on the scale where it 

has the large deviation form 


Pr{x,t) 



( 1 ) 


where positive and negative fluctuations are scaled differently 
and I{y) represents the large deviation function (LDF). For 
systems with x^ —x symmetry, there is only a single length 
scale ^+(f) = ^-(0 = \ {t) and I{y) is symmetric. The LDF 
I{y) characterizes the spatial dependence of the temporal re¬ 
laxation. At different points x, the system relaxes with a 
different x dependent rate. Alternatively, at a hxed point x, 
the system undergoes a dynamical transition from a transient 
to the stationary state at a characteristic time t* (x) such that 
(t*)=x (for positive and negative x respectively). We com¬ 
pute ^±(f) and I{y) explicitly for several systems. We show 
that ^±{t) ^ typically grow algebraically at late times 
with a pair of new exponents. For the underlying stochastic 
process without resetting the typical time dependent length 
scale also grows algebraically ~ where z is the dynam¬ 
ical exponent, e.g., z = 2 for simple diffusion. In presence 
of resetting, we show that v± are generically smaller than z. 
Furthermore, we show that the rate function I(y) exhibits a 
universal feature: it has a pair of singular points at y = 
and y = —yl, corresponding to the frontiers of the core re¬ 
gion —signaling a dynamical phase transition. As discussed 
later, physically, this LDF /(y) and its singularity provides a 
sharp spatio-temporal separation between typical and atypical 
trajectories of the underlying stochastic process with resetting. 
It turns out that, generically, the second derivative /"(y) is dis¬ 
continuous at this pair of singular points. 

Dijfusion of a single particle .— We begin with the simple 
example of a single particle diffusing in one dimension (gen¬ 
eralization to higher dimensions is straightforward) whose po¬ 
sition is stochastically reset to a fixed position (which is taken 
to be the origin without loss of generality) with a constant rate 
r. Let Pr{x,t) be the probability density function (PDF) of the 
position X of the particle at time f. There arise two situations: 
one in which no resetting events occur during the observation 
time f (the probability of which being e^''‘) so that the parti¬ 
cle moves from the origin to x as a free diffusion, and another 
in which the last resetting event before t occurs at time f — T 
(and no resetting occurs in the remaining time T —the PDF of 
which being re^’"^) so that the particle moves as a free diffu¬ 
sion in the hnal stretch of the time T. Therefore, taking into 
account both situations (and integrating over t), one gets ini 


Pr{x,t)=e ''‘PQ{xf)+ [ dire '"^Pq{x,t), (2) 

Jo 


where Poixf) = exp[—x^/ {4-Dt)]/\/4-nDt is the propagator 


for the particle to diffuse freely from the origin to the posi¬ 
tion X in time t in the absence of the resetting. The NESS is 
obtained by taking the t limit in Eq. Q, which yields ||T3| 

Pr{xf^oo) = J^ dxre^’"^Po{x,T) = ^exp{-a\x\), (3) 

where a = \fVjD. To analyse (1^ for hnite f, it is convenient 
to reexpress it using a change ofvariable T = wf, yielding 


Pr[xd) = 


-r<I>(l.x/r) 


ri/f 


dw 


y/AnDt \/AnD Jo 




(4a) 


where 


^{w,y) = rw +. (4b) 

ADw 

For large t and fixed y = xlt, the integral in the second term 
can be estimated by the saddle point method. The function 
<I>(w,y) evidently has a single minimum with respect to w at 
w* = \y\/y/ADr, obtained by setting d„^{w,y)\„^„* — 0. If 
w* < 1, the saddle point occurs within the integration limits 
w £ [0,1] and one gets, from Eq. (|^ Pr{x,t) ~ for 

large f, where 'I>(w*,y) = a |y|. In contrast, for w* > 1, the 
function <I>(w,y) has its lowest value in w G [0,1] at w = 1. 
Hence the integrand in the second term is dominated by the 
regime at w = 1 (and is of the same order as the hrst term). 
Physically, this corresponds to trajectories which have under¬ 
gone zero (or almost zero) resettings up to time t. One then 
gets Pr{xf ) ~ vvith ^(Ijy) = r-|-y^/(4D). Sum¬ 

marizing, we obtain 


Prixf) ^ 




(5a) 


where the LDF 


a|y| for|y|<y*, 

y2 

'■+4D 


(5b) 


with y* = y/ADr. 

Comparing Eq. ( [5al i with Eq. Q shows that the growing 
length scale ^ (t) ~ f, much larger than the typical diffusion 
length scale ~ \/t. The linearity of the LDF for jyj < y* im¬ 
plies that, for any large but hnite f, there is an interior spatial 
region —y*f < x < y*f, where NESS has been achieved, as 
Prixf) ^ exp(—ajxj) becomes independent of t —in agree¬ 
ment with Eq. 0- However, there is still an exterior region 
|x| >y*t that has not been relaxed to the NESS yet. The 
boundaries between the two regions move at a constant speed 
y*. From Eq. (5bi, it is easy to check that while /(y) and 
its hrst derivative are both continuous at y = ±y*, its second 
derivative has a discontinuity at y = ±y*. Therefore, the LDF 
has a second order discontinuity at the points ±y*. Figure]^ 
shows very good agreement between the above large deviation 
form of the PDF and numerical simulations. The above anal¬ 
ysis can be easily generalized to higher dimensional diffusion 
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FIG. 2. (Color online) PDFs of the position of a particle diffusing 
in one dimension with a diffusion constant D = 1 /2, whose position 
is stochastically reset to the origin at a constant rate r = 1 /2. (a) The 
points are simulation data for the PDFs at f = 10 (blue circles) and 
r = 15 (green squares). The (magenta) dashed line is the infinite time 
NESS given in Eq. The vertical dashed lines mark the positions 
±y*?. (b) Same data as in (a) compared with large deviation result 
[Eq. (0] of the PDF (normalized numerically), denoted by the (ma¬ 
genta) dashed lines. The dashed vertical lines at y = ±1 mark the y* 
at which the LDF has a second order discontinuity. 


as well as to other stochastic processes such as the fractional 
Brownian motion (see Appendix). 

What is the physical signihcance of this phase transition? 
The probability density Pr{x,t) can also be interpreted as the 
density at time f of a swarm of independent Brownian mo¬ 
tions, each subjected to stochastic resetting with rate r, all 
starting from the origin at f = 0. Our calculation shows that 
at time t the density for W <y*t becomes stationary, while 
is still time dependent for |x| > y*t. From the analysis above, 
it is clear that, for l-^l > /T the density is typically of the 
form ~ e^''‘Po{x,t) in Eq. (j^, i.e., it corresponds to particles 
that have undergone almost no resetting up to time t. This 
is of course a very rare event and these particles in the outer 
region thus have very atypical trajectories. In contrast, the 
particles in the inner core region correspond to typical tra¬ 
jectories that have undergone a large number of resettings - 
leading to a stationary behavior in this regime. The LDF I(y) 
in Eq. ( |5al l probes precisely the separation between these two 
regions, i.e., between the typical and the atypical trajectories. 
The singularity in the LDF signihes a sharp separation be¬ 
tween these two types of particles. In any typical application 
of resetting, for instance in the optimization of search algo¬ 
rithms, we would ideally like to keep, at any given finite time 
f, only the typical trajectories and not the atypical ones - since 
the latter ones do not feel the resetting at all. The LDF I(y) 
and its associated singularity, that sharply separates the two 
types of trajectories, thus provides a very useful and practical 
way to select the typical ones at any given time f. Even though 
we discuss it here in the context of a single particle diffusion, 
it turns out that this picture is quite generic and holds for ar¬ 
bitrary stochastic processes undergoing resetting and even for 
spatially extended systems, such as fluctuating interfaces that 
we discuss next. 

Fluctuating (1 + 1)-dimensional interfaces .— We next look 
at the effect of resetting on extended correlated objects such 
as a fluctuating interface growing with time over a linear sub¬ 
strate of size L. The interface is characterized by its height 
held H{x,t) which typically evolves via a stochastic dynam¬ 


ics ll22l - [24l . The height fluctuation is measured by the rela¬ 
tive height held, h{xf) — H{x,t) — H{x,t), where H{x,t) = 
^ * fo dxH{x,t) is the spatially averaged height. The rough¬ 
ness of the interface is usually measured by the interface width 
W{L,t) defined as W^{L,t) = ^ J^dx (h^(x,t)). For a large 
class of interfaces, W{L,t) increases with time as W(L, f) ^ f ^ 
for 0 ^ f <C L^, before saturating to a time-independent value 
Wsat ~ for t ^ U-. The scaling exponents a, p and z are 
known as the roughness exponent, growth exponent and dy¬ 
namic exponent respectively, and are related by the scaling re¬ 
lation z = a/P, leaving only two independent exponents 
Moreover, in the growing regime 0 <^t <^U-, the full height 
distribution Po{li,t) has the generic scaling form ll22l - E4l 

Po{h,t)z^{Tt)-Pg{{Tt)-Ph) , ( 6 ) 

where the scaling function g{x) is identical for all models be¬ 
longing to the same universality class while F^' is a model- 
dependent microscopic time scale. For example, for simple 
linear stochastic interface models belonging to the Edwards- 
Wilkinson (EW) universality class (where P = 1/4 and z = 2), 
the height distribution at all times is simply Gaussian ll25ll 
with g(x) exp[—x^]. Another widely studied class of grow¬ 
ing interfaces in (1 -F 1) dimensions belong to the Kardar- 
Parisi-Zhang (KPZ) universality class Il22lj24l |2^ . In this 
case one has p = 1/3 and z = 3/2 while the scaling func¬ 
tion g{x), for a flat initial condition H{x,t = 0) = 0, is related 
to the Tracy-Widom distribution ll27l associated to the Gaus¬ 
sian Orthogonal Ensemble (GOE) of random matrices ||28l, 
/i(z)’ which describes the fluctuations of the largest eigen¬ 
value of GOE matrices. One has indeed g{x) = /i(x+ (x)) 
where (x) = J//’ocdXfiix)- While the full form of fi{x) 
is rather nontrivial— the tails have simpler non-Gaussian 
forms: fi{x) ^ exp(-|xp/24) as j ^ -oo and fi{x) 

exp(-2x3/2/3) as j oo. 

Let us now consider the height held of such a generic (1 + 
1) dimensional interface evolving under its own dynamics and 
subject it to resetting to its initial height profile at constant rate 
r. Following our discussion prior to Eq. Q, one can relate 
the height distribution Pr{h,t) in the presence of resetting to 
that of Pq {h,t) without resetting via the same equation Eq. (0, 
with X replaced by h. Using the scaling form for Po{li,t) in 
Eq. (0 and making the change of variable T = wt, Eq. Q 
reduces to 

Pr{h,t)^{rt)-Pe-’-‘g{{rt)-Ph) 

4- rt{rt)^^ J (7) 

As before, Pr{h,t) approaches a stationary distribution as f —>■ 
oo for any r > 0 EOl 

Pr{h,t ^ « iT/r)-I^Gp {{T/r)-Ph), (8a) 

where the scaling function is given by 

Gp {x) = ^ dyy-f^e-^'g{xy-P). 


(8b) 
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FIG. 3. (Color online), (a) The PDFs of the relative heights of a 
periodic interface of size L = 2*^ at various times, generated from a 
TASEP (shown by points) are collapsed by choosing F = 2.4 to the 
TW GOE PDE with zero mean, shown by the (magenta) dashed line, 
(b) The points are simulation data for the PDF of the relative heights 
at f = 50 (blue circles) and t = 100 (green squares), with resetting 
with rate r = 0.05. The (magenta) dashed line is the infinite time 
NESS given in Eq. l[^. The vertical dashed lines mark the positions 
(c) Same data as in (b) compared with Eq. 1 16 1 
denoted by the (magenta) dashed lines. The vertical dashed lines 
mark the positions • (4) The LDFs computed numerically from 
jl6|l for various t are compared with (|19b|>. 


Now to investigate the approach to the NESS, we consider 
the generic case when g{x) ^ exp(—a±|x|^±) as x —)■ ±oo. For 
example, for the KPZ with flat initial condition, 7 + = 3/2, 
a+ = 2/3 and 7 _ =3, = 1/24. Substituting these generic 

tails of g(x) in Eq. ( [T^ we obtain that the large deviation form 
of Pr{h,t) (see Appendix) 


Pr{h,t) 




where the LDF is given by 


r 

I{y) = } pv±{yiy^ 

I r + (7±br* 


for b| < 7 i, 
for b| > yi . 


(9a) 


(9b) 


The exponents v± = 7±/(l + j37±) and the boundaries y^z 
as well as the constants b± can be explicitly computed (see 
Appendix). The signs ± are chosen for y > 0 and y < 0 re¬ 
spectively. Comparing Eq. (19ai and Eq. ([T]), we see that the 
growing scale, separating the inner NESS regime from the 
outer transient regime in the height space, ^±{t) ~ is 
asymmetric (different respectively for positive and negative 
height fluctuations). The height fluctuations in the intermedi¬ 
ate range —y* f < h < y^f reaches NESS. Moreover, 
/(y) again has a second order discontinuity at the two singular 
points y%. 

To perform simulations, we consider an interface model 
generated from the well known totally asymmetric simple ex¬ 


clusion process (TASEP) 130) on a ring. This interface model 
is known to belong to the KPZ class ll22l l24l . Figure (a) 
shows that the scaled PDF for various times in the growing 
regime (1 ^ f ^ U-) can be collapsed by choosing an unique 
value of r (same for all f), to the TW GOE PDF with the mean 
shifted to zero. We use this value of E subsequently. We next 
simulate the dynamics of the interface in presence of resetting 
(with rate r) to the flat initial condition. It is evident from 
Fig. [^(b) that at long times the central part \h\ <C of the 
distribution of the relative heights reaches the NESS given by 
Eq. ([T^. Figure|3(c) shows very good agreement between the 
simulation and the theoretical results obtained from numerical 
integration of Eq. In Fig. 1^ (d) we plot the LDFs com¬ 
puted numerically from ([T^ for various f, which converges to 
( [T9bl i as time increases. 

Discussion. — In summary, we have studied analytically the 
approach to the stationary state of several systems subjected 
to resetting to the initial condition at a constant rate. We 
have shown that the relaxation mechanism in these systems 
is highly unusual. In typical systems approaching to their sta¬ 
tionary states, the late time relaxation is governed by a single 
time scale independent of space. In systems with resetting 
studied here, the late time relaxation is space dependent and 
is characterized by the LDF I(x/^ (t)) defined in Eq. Q. The 
growing length (or height) scale around the resetting center 
behaves at late times as ^ (f) ^ f where v < z is a new ex¬ 
ponent typically smaller than the dynamical exponent z of the 
process without resetting. For system without the x —> —x 
symmetry, such as in the KPZ equation, there are actually a 
pair of asymmetric length scales ^±(f) ^ ?*/''=•= where v± < z. 
We have computed this LDF I{j) explicitly in several exam¬ 
ples and found that it exhibits a universal feature; it has a 
a pair of singular points where the second derivative is dis¬ 
continuous. The singularities in I{y) provide a sharp separa¬ 
tion between typical and atypical trajectories and thus can be 
useful in various applications such as search optimization us¬ 
ing resetting, in detecting and getting rid of the undesirable 
atypical trajectories (that have not undergone resetting) at any 
finite time t. The important point is the presence of these sin¬ 
gularities providing a sharp separation - the actual order of 
the singularities is of less relevance. We conclude by noting 
that singularities in the LDFs have also been found recently 
in several different contexts and there is a growing interest in 
understanding the significance of these singularities liiHiTi. 
In this paper we have provided a clear physical meaning of 
these singularities that act like separatrices between typical 
and atypical trajectories. 

Acknowledgment. — We acknowledge support of the Indo- 
French Centre for the Promotion of Advanced Research under 
Project No. 4604-3. 

Appendix. — In this appendix, we give some details of the 
calculations described in the main text. 

Fractional Brownian Motion. — The analysis of Pr{x, t) car¬ 
ried out in the main text above for the simple Brownian mo¬ 
tion can be easily generalized to the fraction Brownian motion 
(fBM), which represents a Gaussian process with zero mean 
and two-time correlator, (x(fi)x(f 2 ) = t\^ +^ 2 ^ ~ \h 
where 0 <H < \ represent the Hurst index. Clearly FI =1/2 
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corresponds to the simple Brownian motion. For simplicity 
we discuss the one dimensional case, though generalization to 
higher dimensions is straightforward. The propagator for the 
fBM in one dimension, starting initially at the origin, again 
has the simple Gaussian form 


PQ{x,t) 


1 

, exp 


4Dt^^ 


( 10 ) 


We then subject the particle undergoing fBM to resetting at the 
origin with a constant rate r. Then, the propagator Pr{x,t) in 
presence of the resetting satisfies the generic relation given by 
Eq. (|^. Substituting the bare propagator Po{x, t) from Eq. ( [loj i 
in Eq. Q, we can then carry out the same analysis as the H ~ 
1 /2 case in the main manuscript. As in the case of // = 1 /2, it 
turns out that for general 0 < H < 1, the growing length scale 
^ (f) ^ much larger than the typical spread ~ for 

large t. In the scaling limit x ^ t ^ oo with fixed, 

we obtain, using the saddle point method (discussed already 
for H = 1/2 in the main manuscript) 


where the large deviation function (LDE) is 

for|y|</, 

/ \ \ f H \ ^ 

witharr = r(^l + -j(^—j and/= 



(11a) 


(lib) 



It is easy to check that both the LDE and its first derivative are 
continuous across ±y*, while the second derivative has a dis¬ 
continuity at ±y*. The NESS is realized in an interior region 

where Pr{x,t) ^ exp(—|jic/+'/2) becomes independent of 
time, with the boundaries between this NESS and the outer 
transient regime moving out with time as \x* \ = 

The case of the Kardar-Parisi-Zhang equation .— A widely 
studied class of growing interfaces in (1 -f 1) dimensions be¬ 
long to the Kardar-Parisi-Zhang (KPZ) universality class ll22l - 
l24l . Here the height field evolves via the nonlinear KPZ equa¬ 
tion ll26l 


dH d^H 
dt (9/ 


dH 

dx 


rj{x,t), 


( 12 ) 


where v is the surface tension, A represents the strength of the 
nonlinearity, and rj{x,t) is a Gaussian white noise with zero 
mean and correlations {rj{x,t)rj{x\f)) = 2D5{x — x')5(t — 
t'). In absence of nonlinearity (A = 0), the KPZ equation (12 1 
reduces to the EW equation. 

The scaling exponents for the KPZ equation in one dimen¬ 
sion are well known Il22l : a = 1/2, j3 = 1/3 and z = 3/ 2. 
In the KPZ case, the spatially averaged height H{xf) = 
P ^ fo dxH{x,t) grows linearly with time with a non-zero ve¬ 
locity Voc = {Z/2)L^'- Jq {{dH/dx)^) dx. In contrast, for the 


EW case, H{x,t) ~ y/tJZ for large t. In the growing regime 
0 ^ f ^ to which we focus below, the PDE Po{h,t), while 
trivially Gaussian for the EW case, is highly nontrivial in the 
KPZ case and has been a subject of intense investigations in 
recent times. It has been solved exactly for the KPZ equation 
only very recently ll28l . It turns out to depend on the initial 
condition of the height profile 12811^ . 

Eor example, for aflat initial profile, the height H{x,t) of 
the interface can be written as 

H{x,t)=v^t + {Tty/^X{x), (13) 

where E is a constant that depends on the parameters of 
the interface model, and is a time independent random 
variable distributed according to the so-called Tracy-Widom 
(TW) distribution fi{x), that characterizes the fluctuations 
of the largest eigenvalue of random matrices in the Gaus¬ 
sian orthogonal ensemble (GOE) l27l . While the full form 
of fi{x) is rather nontrivial— the tails have simpler non- 
Gaussian forms: f\{x) ~ sxp(—|;|^p/24) as x ^ and 
/i (z) 6tcp(—2x^/^/3) as X —>■ oo. In terms of j, we have 

h{x,t) = {Tty^^[x{x)-X], (14) 

where Z = A * /o' X{x) dx. The law of large number dictates 
that X {%) lit the limit °°,so that {h) = 0. Therefore, in 
this case, in the scaling limit of t ^ h ^ while keeping 
hjfll'^ fixed, the height fluctuation is distributed according to 

p^{hf) «(rf)-'/3/i {{Tt )-'/\), (15) 

where fi{x) = fAx+ (z))- This is indeed of the form an¬ 
nounced in Eq. ® in the main text above with j3 = 1/3 and 
g{x) = f\{x). Similarly, for the droplet initial configuration 
(curved geometry), the scaling function for the height distri¬ 
bution is the TW distribution corresponding to the Gaussian 
unitary ensemble (GUE). 

Generic interfaces: saddle point calculation. — Let us now 
analyze the height held of a generic (1 -f 1) dimensional inter¬ 
face evolving under its own dynamics and subject it to reset¬ 
ting to its initial height profile at constant rate r. In this case, 
the full height distribution P(){h,t), in the absence of resetting 
(r = 0), has the generic scaling form ll22l - [24l given by Eq. (j^. 
As shown in the main text, one can relate the height distri¬ 
bution Pr{h,t) in the presence of resetting to that of Po{hf) 
without resetting via the Eq. Q with x replaced by h. Using 
the scaling form for Po{h,t) in (j^ and making the change of 
variable T = wt, it reduces to 

+rt{rt)^l^ J (16) 

In this case, as in the case of single particle diffusion, Pflhf) 
approaches a stationary distribution as t ^ for any r > 0, 
given by ll^ 

Pr{h,t ^ « (r/r)-%p {{P/r)-^h), (17a) 
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where the scaling function is given by 

Gj 3 (x) = ^ dyy^^e^-g{xy^^). (17b) 

Now to investigate the approach to the NESS, we consider 
the generic case when g(x) ^ exp(—a±|x|^±) as ±oo. For 

example, for the KPZ with flat initial condition, 7 + = 3/2, 
= 2/3 and 7 _ = 3, a_ = 1/24. Substituting these generic 
tails of g(x) in Eq. ( [T^ we obtain 

Pr(h,t) - 

+ rf(rf)^^ (18a) 

Jo 

where v± = 7±/(1 + j37±) and 

, N b±\y\'^^ <^± 

4>(w,v) = rwH-s- with b± = —j ,—. (18b) 


With respect to the variable w, the function <I>(w,y) has 
a unique minimum at w* = (jyl/y^) where y!/ = 

(r/[j 37 ±Z 7 ±])'^^*. As in the case of single particle diffusion 
discussed after Eq. (4bi, if w* < 1, the most dominant con¬ 
tribution to Eq. (18a I comes from the neighborhood of w* so 
that Pr{h,t) ~ whereas for w* > 1 , it is domi¬ 

nated by the boundary terms so that Pr{h,t) ^ 

Therefore, the large deviation form of Pr{h,t) is given by 


Pr{h,t) 




where the LDF is given by 


ny) 


r\y\^^ 

j3v±(yi)v± 

r + b± \y\'^^ 


for |y| <yi, 
for |y| >yi, 


(19a) 


(19b) 


where the ± signs are chosen for y > 0 and y < 0 respectively. 
This yields the expression given in Eq. (9b), together with 
the explicit expression for the boundaries yj. and for the con¬ 
stants b± [see Eq. ( |18b 1 and below]. Moreover, one can show 
that /(y) in ( 19b| has a second order discontinuity at the two 
singular points y^. 

Generic second order discontinuity in the large deviation 
function .— We now consider a generic system described by 
a stochastic variable x (which may represent the position of a 
particle undergoing generic stochastic dynamics, or may rep¬ 
resent the height of a fluctuating interface). Let Po{x,t) be the 
bare probability density function (PDF) of x at time f, starting 
at .r = 0, in the absence of resetting. Let Pr{x,t) denote the 
PDF of X at time t in the presence of resetting to the initial 
condition .r = 0 at a constant rate r. The PDF Pr{x,t) is re¬ 
lated to the bare propagator Po{x, t) via the general relation in 
Eq. (j^. In the limit t ^ Pr{x,t) in Eq. ([^ approches the 
stationary distribution 


Pr(x,t^o°) = r f die '^^Po{x,t), 
Jo 


( 20 ) 


which is fully determined by the bare PDF Po{x,t). 

In order to study the approach to the stationary state at late 
times, we need to analyze Eq. 0 for finite but large t. To 
proceed, it is convenient to make a change of variable i = wt 
in the integral which yields 

Pr{x,t) = e^'''PQ{xf) + rt [ dwe^*™ Po{x,wt). (21) 

Jo 

We are interested in the behavior of Pr{x,t) for large x and 
large t. Hence, inside the integral we need to substitute the 
large x, large t behavior (with w fixed) of PQ{x,wt). For 
generic self-affine systems where x{t) ^ f^, the PDF Po{xf) 
is expected to have a scaling form, for large x, large f, keeping 
x/tP fixed 


Po{x,t) : 


irt)P^\(rt)P 


( 22 ) 


where g{y) is the scaling function and E^* is a model- 
dependent microscopic time scale, not important for our anal¬ 
ysis, and henceforth set to unity without any loss of generality. 
Substituting the scaling form Eq. (22i of Po^xf) in Eq. •HD 
gives 


Pr{x,t)=t Pe 


-rt 
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tp 
1 dw 
wP 


e-‘’-'^g 


\{wt)P 


(23) 


Next, we anticipate that Pr(xf) will have a large deviation 
form [see Eq. Q] 

P,(x,f) - exp (-tl > (24) 

where /(y) is the rate function and ^ (f) represents the growing 
length scale associated with atypically large fluctuations of x 
that are much bigger than the typical fluctuations x^tP. Note 
that in some cases (such as for the fluctuating interfaces be¬ 
longing to the KPZ universality class), the x —> —x symmetry 
is broken. In such cases, the positive and negative large fluc¬ 
tuations of X occur at different length scales it) and accord¬ 
ingly, Pr{x,t) — exp[—f/(.r/^±(f))] for positive and negative 
fluctuations respectively [see Eq. Q]- In these systems, we 
then need to scale the positive and negative fluctuations dif¬ 
ferently. However, in problem with x^ —x symmetry, there 
is only a single scale i*(f) and the LDF /(y) is symmetric. 
For simplicity, we restrict ourselves below only to the posi¬ 
tive fluctuations, as the analysis of the negative fluctuations is 
similar. 

In order to probe the behavior of Pr{x,t) on the scale x — 
^ (f) ^ f^, we then need to substitute the tail behavior of g{y) 
for large y in the integrand of the second term in Eq. (j2^. 
Generically, g{y) — exp[—ay^ for large positive y, where a is 
an unimportant model dependent constant. Substituting this 
tail behavior in the integrand in Eq. (j2^ we get 

+ rt^-P dww-Pe-‘^^'^’^^^P^\ (25) 

Jo 
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FIG. 4. Schematic plot of the function ^{w,y) as a function of 
w (but fixed y), which has a unique minimum at w*{y), determined 
from Eq. ([27j for the two cases (a) w* < 1 and (b) w* > 1. For the 
case (a), the integral in Eq. l |25[ l is dominated by the minimum w*, 
whereas for the case (b) the integral is dominated by the upper limit 
w = 1 of the integral. The minimum values of <F(w,y) in the interval 
w e [0,1] in both cases, are marked by the red circles. 


where ^ (f) = ^ f/3 

<I>(w,y) = rw- 


ay' 


(26) 


Evidently as a function of w (but fixed y), has a 

unique minimum at w*(y), determined from 




= 0 where = 


dw 


(27) 


Actually the subsequent analysis will be very general and we 
do not need to use the specific form of <I>(w,y) —the only fact 
we will use is that <I>(w,y) has a unique minimum at w*{y). 


determined via the minimization in Eq. ( |27| ). 
Eor large t, the integral in Eq. 


is dominated by the min¬ 
imum w*(y), as long as w*{y) < 1, i.e., the saddle point lies 
within the limits of the integral w S [0,1], as shown in Fig.|^ 
For w*(y) > 1, the integral is dominated by the upper limit 
w = 1 of the integral (see Fig. |^, which is of the same order 
as the first term. Thus, we have a critical point at y = y* which 
is determined from w*(y*) = 1. Thus evaluating the integral 
in Eq. ( p5| ) by the saddle point method for large t gives the 
large deviation form given in Eq. (p4jl where 


iiy) = 


^(w*(y),y) fory<y*, 


^ihy) 


for y > y* . 


(28) 


This is the generic mechanism of the dynamical transition— 
the rate function /(y) changes its behavior at y = y*. While 
the function /(y) is evidently continuous at y = y*, we show 
below that while its first derivative /'(y) is also continuous at 
y = y*, the second derivative I"{y) is discontinuous, signaling 
a second order dynamical transition. 

Consider first a function of y of the form H{y) = <I>(M(y) ,y) 
where the first argument M(y) depends implicitly on y. The 
derivatives of H{y) can be easily determined by the chain rule. 
For example, the first derivative is given by 

H'iy) = ^ = du‘i>{u,y) u'{y) + dy^{u,y). (29) 


Similarly, the second derivative is given by 

H"(y) = [u'{y)Y + 2dudy^{u,y)u'{y) 

+ <9„<!>(«,y) u”(y) + dyMu,y) ■ 


(30) 


Now, let us first consider the case y < y* in Eq. (28 1 . To 
evaluate I'{y), we substitute u{y) = w*(y) in Eq. (29i. Using 
the minimization condition in Eq. ( (27) i, we obtain 

/'(y) =ay(I)(w*(y),y) for y < y*. (31) 

If now y —>■ y* from below, using w* (y*) = 1, we get 

= (32) 

.y=r 


Fory >/ from above, we have /(y) = <I>(l,y) from Eq. (Cm. 

Hence as y —>■ y* from above, I'{y ^ y*) = <9j.O(l,y) 

_ v=y’^ 

Comparing this with ( |32] l, we see that /'(y) is continuous at 

% 

y = y . 

We next consider the second derivative I" (y) as y —y y* from 
below and above. Consider first the case y </■ Substituting 
u{y) = w*(y) in ( |30l l and using the minimization condition 
Eq. we get 

7"(y) = ,y) [(w*)'(y)]^ 

+ 2d„*dy^(w*,y) (w*)'(y) + 0(w*,y). (33) 

Now, as y —>■ y* from below, w* (y) —> 1 and we get 

7"iy^y*)= dl,^{w*,y) [(w*)'(y)]^ 

+ 2d„*dy<P{w*,y) (w*)'(y) 

+ dMhy) ■ (34) 


VP* = 1,);=)?* 


In contrast, for y > y*, /(y) = <I>(l,y) and hence I"{y) = 
dy<i>{l,y). Hence, as y —>• y* from above, /"(y —y y*) = 

(9,^<I>(l,y) . Comparing this with Eq. (34i, we see that 

>•=>'* 


generically the second derivative of the ratemnction is dis¬ 
continuous across y = y* and the value of the discontinuity is 
given by 


lim [/"(y*-e)-/"(y* + £)] = [dlM^\y) [(w*)'(y)]' 


+ 2d„.dy^{w\y){w*)'{y) 


w* = l,3?=y* 


(35) 


As an example, let us consider the case of diffusion with 
resetting (the first example in the main text). In this case, from 
Eq. (|g 

,.2 


<I>(w*,y) = rw* -F 


T 


ADw* 


(36) 


where w*(y) = yjy/ADr (considering only the positive side). 
It is then straightforward to evaluate the discontinuity in 
Eq. ( p5| ) and we get 

lim [I"{y* - e) - I'\y* + e)] = - A . (37) 


2D 





























Similarly, one can obtain the value of the discontinuity in the second derivative for the other examples discussed in the main 

text. 
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